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We study the stability of the pentaquarks uuddQ, uudsQ and udssQ (Q 



I/-} ' c, b, or t) of positive parity in a constituent quark model based on Goldstone 

boson exchange interaction between quarks. The pentaquark parity is the 



antiquark parity times that of a quark excited to a p-shell. We show that the 



> 

■^J- ■ Goldstone boson exchange interaction favors these pentaquarks much more 

m ' 

O 

OO 



than the negative parity ones of the same flavour content but all quarks in 
the ground state. We find that the nonstrange pentaquarks are stable against 
strong decays. 

The existence of particles made of more than three quarks is an important issue of QCD 
inspired model. The interest has been raised so far by particles described by the colour state 
[222] the tetraquarks q 2 q 2 , the pentaquarks q 4 q and the hexaquarks q 6 . The present study 
is devoted to pentaquarks, first proposed independently by Gignoux, Silvestre-Brac and 
Richard Jl| and Lipkin about ten years ago. Within a constituent quark model based on 
one-gluon exchange (OGE) interaction, these authors found that the states P® s = \uuds~c) 
and = \uddsc) and their conjugates are stable against strong decays. Within better 
approximations, they turned out to be unstable ||f|]. A systematic theoretical study in 
a model with OGE interaction suggested several candidates for stability, and among them 
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especially those with strangeness S = -1 or -2. In particular, the uudsc system was bound 
by -52 MeV. The nonstrange systems uuddQ (Q = c or b) were unbound. 

Calculations done within an instanton model || or a Skyrme model [[/j] show that pen- 
taquarks, irrespective of their strangeness, appear as bound or near-threshold resonances 
depending on the model parameters. Moreover, the lowest pentaquark states predicted in 
Ref. have positive parity (L = 1). 

If bound, the lifetime of the pentaquark uudsc or uddsc is expected to be similar to 
that of the Df meson. Using various mechanisms, the typically estimated pentaquark 
production cross section is of the order of 1 % that of Df || . Based on the above theoretical 
predictions, experiments are being planned and the first search for the pentaquarks and 
P z ~, performed at Fermilab, has just been reported ||. The decay P° s — > (pnp was analyzed 
for two different pentaquark masses : M = 2.75 GeV, being the lowest mass expected from 
the OGE model and M = 2.86 GeV, the value at which the largest number of events was 
observed. No convincing evidence for pentaquarks decaying to (pirp was observed yet. 

The theoretical predictions are definitely model-dependent. In a previous work | |10| . 
we studied the stability of heavy-flavoured pentaquarks within a chiral constituent quark 
model |TT|-fn|] originally proposed by Glozman and Riska [|ITJ. In this model, the hyperfine 



splitting in hadrons is obtained from the short-range part of the Goldstone boson exchange 
(GBE) interaction between quarks, instead of the OGE interaction of conventional models, 
as discussed above. The main merit of the GBE interaction is that it reproduces the correct 



ordering of positive and negative parity states in all parts of the considered spectrum [|T^[T4 
in contrast to any other OGE model. 



In Ref. ||10|| , we considered pentaquarks with strangeness ranging from S = -3 to S = 
0. There, we assumed that all light quarks are identical and the ground-state orbital (O) 
wave function is symmetric under permutation of light quarks, i.e. it corresponds to the 
Young diagram [4] . The subsystem of light quarks must necessarily be in a colour 3, or 
alternatively a [211] c state. Then, the Pauli principle allows a certain number of spin [f] s 
and flavour [f] F states to be combined with [4] Q and [211] c to give a total antisymmetric 



four-quark state. We found that any of these states combined with a heavy antiquark Q 
where Q = c or b gave rise to a pentaquark energy which was at least 300-400 MeV above 
the dissociation threshold nucleon plus meson, i.e. the considered pentaquark cannot be a 
bound compact object. Its parity is P = -1, due to the antiquark. 

The novelty of this study is that, within the same GBE model, we analyse the stability 
of pentaquarks with P = +1. In such a case, the parity of the pentaquark is given by 
P = (—)"'" \ thus the light quarks must carry an angular momentum L odd. Here, we 
consider the case L = 1, which implies that the subsystem of four light quarks must be in a 
state of orbital symmetry [31] G . Although the kinetic energy of such a state is higher than 



that of the totally symmetric [4] Q state, a schematic estimate suggests that the [31] G 
symmetry would lead to a stable pentaquark. In the following, we give the arguments of 



Ref. |14| based on a simplified GBE interaction of the form 



C X Yj X i - X j ff i' ff 3 



(1) 



i < j 



with Aj the Gell-Mann matrices, B{ the Pauli matrices and C x = 30 MeV, determined from 



the A-N splitting [11 



In the spirit of Glozman and Riska's model, there is no meson-exchange interaction 
between quarks and antiquarks. It is assumed that the qq pseudoscalar pair interaction is 



automatically included in the GBE interaction |TT| . Then, as far as the spin-spin interaction 
is concerned, the discussion is restricted to the light q 4 subsystem. The Pauli principle 
allows for the following two totally antisymmetric states with [31] symmetry, written in 
the flavour-spin (FS) coupling scheme 



|1> = ([31] [211] c 



oc 



[22] F [22] 5 [4] FS ) 



(2) 



|2> 



[31] [211] 



c 



oc 



[31] F [31] 5 [4] 



FS 



(3) 



Asymptotically, a ground state baryon and a meson, into which a pentaquark can split, 
would give [3] x [2} = [5] Q + [41] + [32] Q . By removing the antiquark, one can make 



the reduction [41] — > [31] x [1] G or [32] — > [31] x [1] Q . Thus, the symmetry [31] 
of the light quark subsystem is compatible with an L = 1 asymptotically separated baryon 
plus meson system. 

The expectation value of (1) calculated, for example, according to the Appendix of Ref. 



|T5|| , is —28 C x for |1) and —64/3 C x for |2). These two states would actually couple via 



1 3 

a quark-antiquark spin-spin interaction to a total angular momentum J= — or — where 
J = L + S + sq, with L, S the angular momentum and spin of the light system 
and sq the spin of the antiquark. As the quark-antiquark interaction is neglected here, in 
the following we restrict our discussion to the lowest state, i.e. 1 1) . The quark-antiquark 
interaction is neglected in the description of mesons as well, as for example in Ref. ||16|| , so 
that the meson Hamiltonian contains a kinetic and a confinement term only. 
We are interested in the quantity 

AE = E(q A Q) - E(q 3 ) - E(qQ) (4) 

In our schematic estimate, we suppose that the confinement energy roughly cancels out in 
AE. Then, the kinetic energy contribution AKE to AE is AKE = 5/4 hu in a 
harmonic oscillator model and the GBE contribution for the state 1 1) is AV X = —14 C x . 
With Tiu m 250 MeV, determined from the N(1440) - N splitting JTTJ] , this would give 
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AE = - hu - 14 C x = - 107.5 MeV (5) 

i.e. a substantial binding. This is to be contrasted with the negative parity pentaquarks 
studied in Ref. |L(| where for the lowest state one has AE = 3/4 Uuj — 2 C x = 127.5 MeV, 



i.e. unstability, consistent with the detailed study made in [fL0| . 

The estimate (5) is a consequence of the flavour dependence of the GBE interaction. For 
a specific spin state [f] s , a schematic OGE interaction of type V c m = — C c m A?.A^ <?j. (re- 
does not make a distinction between [4] G and [31] so that the [31] state will appear higher 
than [4] G due to the kinetic energy. The GBE interaction overcomes the excess of kinetic 
contribution in [31] G and generates a lower expectation value for [31] than for [4] . 



The GBE Hamiltonian has the form |12P : 



H = Em + j:£r- %f^r + E v^m) + E vAm) , (6) 

i i Z '" L i *l^i"H j<j iK j 

with the linear confining interaction : 

V ooot (r ij )=~^-X C j Cr ij , (7) 

and the spin-spin component of the GBE interaction in its SUf(3) form : 

V x( r a) = IE K(ry)Af Af 
If=i 

+ E ^(^)Af Af + ^(r y )A?A5 + K/fc^A?) * • a„ (8) 



F=4 



with A = y 2/3 1, where 1 is the 3x3 unit matrix. The interaction (8) contains 7 = it, K, r] 
and rj' meson-exchange terms and the form of (r^) is given as the sum of two distinct 
contributions : a Yukawa-type potential containing the mass of the exchanged meson and a 
short-range contribution of opposite sign, the role of which is crucial in baryon spectroscopy. 
For a given meson 7 = n, K, rj or rj', the meson exchange potential is 



V,(r) = ^T^—iOir - r Q )tf— - ^ exp(-a 2 (r - r ) 2 )} (9) 
For the Hamiltonian (6)- (9), we use the parameters of Refs. fl2| , |T3"|| . These are : 

9jrq_ _ 9rjq_ _ 9Kq _ g gy 9rfq _ ^ 

An Air Air ' An 

r = 0.43 fm, a = 2.91 fm' 1 , C = 0.474 fm~ 2 , m u4 = 340 MeV, m s = 440 MeV 

ji n = 139 MeV, jL n = 547 MeV, ji v , = 958 MeV, fi K = 495 MeV. (10) 



The masses of the threshold hadrons are calculated variationally as in Ref. |10| where 
we assume an s 3 configuration for baryons. They are given in Table 1 where a theoretical 
ts meson is also included in order to analyse the large rriQ limit. The mass of the ts meson 
has been determined by taking m t = 175 GeV. 

For pentaquarks, we used the internal Jacobi coordinates 



x = fx — f 2 , y = (fi + r 2 - 2f 3 ) / 

(11) 

<? = (ri + r 2 + r 3 - 3r 4 ) /VQ , t = (ri + r 2 + r 3 + r 4 — 4f 5 ) / \/lO 

First, we expressed the g 4 orbital wave functions of symmetry [31] in terms of the above 
Jacobi coordinates. We assumed an s 3 p structure for [31] G and inspired by Moshinski's 
method PT7f] , we found the content of the three independent [31] G states flBfl , denoted below 
by ipi, in terms of shell model functions \n £ m). The result is 
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(x|000) (y |000) (z|001) 



(x |000) (y |010) (z|000) 



(f |010) (y |000) (f |000) 



(12) 



(13) 



(14) 



where, for convenience, we took the quantum number m = everywhere. The pentaquark 
orbital wave functions ipf are then obtained by multiplying each ipi by the wave function 
^t|000^ which describes the relative motion of the q 4 subsystem and the antiquark Q. 
Assuming two variational parameters, a for the internal motion of q 4 and b for the relative 
motion of q 4 and Q, we have explicitly 



ipf = N exp 



a ( 2 2 2\ " .2 

~ [ x + y + z ) - ?t 



z y in (z 



(15) 



ipl = N exp 



x z + y z + 2 Z - - t 



yYi 
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(16) 



ipl = N exp 



- - (x 2 + y 2 + z' 



x Y" 10 (a 



(17) 



where 
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N 



2 3/2 a ll/4 6 3/4 



(18) 



The kinetic energy part of (6) can be calculated analytically. For the state (2) or (3), its 
form is 

1 



(T) 



(19) 



ft2 ( 11 
n a 

I 2^ 



with 



4 

where gi, g2 are light quarks and 



2/i 2 



1 3 3 

1 /or gig 2 

mi m2 



(20) 



? 2 

i ^ 2 2 

— + — /or 9i9 2 

I, mi m 2 



5 1 4 
/i 2 /ii m Q 



(21) 



mg representing the heavy antiquark mass. Here, we choose m c = 1.35 GeV, = 4.66 



GeV according to Ref. |T(| and m t = 175 GeV. Taking m u = = m s = tuq and setting 
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a=b, the identical particle limit (T) = - huj with hcu = 2 ah jm is recovered correctly. 



Integrating in the coulour space as shown in Ref. [10] , the confinement part of (6) becomes 

C 



(Vconf) 



(6 (r 12 ) + 4 (r 45 )) 



(22) 



where the coefficients 6 and 4 account for the number of light-light and light-heavy pairs, 
respectively, and 



(23) 



An analytic evaluation gives 



20/1 
9 V Tra 



(24) 



and 

^ = 375 [ 2 \R + ^ (I + 1 ) (25) 

The matrix elements of the spin-flavour operators of (8) have been calculated using the 
fractional parentage technique described in Ref. [113] based on Clebsch-Gordan coefficients of 



the group £4 ||19|| . In this way, each matrix element reduces to a linear combination of two- 
body matrix elements of either symmetric or antisymmetric states for which Eqs. (3.3) of 
Ref. dy] can be used to integrate in the spin-flavour space. The resulting linear combinations 



contain orbital two-body matrix elements of type (ss \V 7 \ ss) , (sp \V^\ sp) and (sp |V^|ps) 
where 7 = tv, K, 77 and rj as in Eq. (8). 

In Table 2, we present results for S = 0, S = -1 and S = -2 pentaquarks. The quantity AE, 
defined by (4) and exhibited in column 5, is obtained from E (q 4 Q^j calculated variationally 
for the state |1) defined by (2). The optimal values of the parameters a and b are indicated 
in each case. In all cases, one has a > b. The inverses 1/a and 1/b give an idea of the 
quark-quark and quark- ant iquark distances, respectively. One can see that at equilibrium 
the light quarks are clustered together, orbitting around the heavy antiquark. 

The present variational solution does not give binding for uudsQ and udssQ. The quan- 
tity AE decreases smoothly with increasing tuq but does not become negative even for 
m t =175 GeV. However, as expected from the discussion around Eq. (5), there is much less 
"repulsive" effect in these type of pentaquarks than that appearing for the negative parity 
pentaquarks studied in Ref. | I~0|| , where AE ~ 370-490 MeV for pentaquarks with similar 



flavour content. 

But the nonstrange positive parity pentaquarks uuddQ are bound by -75.6 MeV, -95.6 
MeV and -102.8 MeV for Q = c, b and t, respectively. The reason is that the GBE interaction 
is stronger in the nonstrange case because of the l/miUij factor in Eq. (9). Thus, the GBE 
model suggests that the nonstrange positive parity pentaquarks are the best candidates 
for stable compact systems. This in contrast to OGE based models where strangeness is 
required || in order to reach stability for heavy-flavoured pentaquarks. The present results 



have similarities with those obtained in J?j from the Skyrme model : 1) the lowest pentaquark 
states have positive parity for any flavour content ; 2) stability does not necessarily require 
strangeness. 

From Table 2 we can obtain upper limits for the masses of the stable quarks as M(uuddc) 
= 2.895 GeV and M{uuddb) = 6.176 GeV. 

It would be very useful indeed to look for a better variational solution for positive parity 
pentaquarks to see whether the soft repulsion persists or disappears for strange systems and 
to improve the upper bounds for the masses of the stable states. A possible way would 
be to consider interference effects with other asymptotic channels, as for example, in Ref. 
@. It would be also interesting to see how the binding energy of the pentaquarks depends 
on the particular parametrization of the hyperfme interaction (nonrelativistic |T2|JT3|| versus 
semi-relativistic GBE model). 



Irrespective of the type of hyperfme interaction, an important issue remains the model 
of confinement. Both in OGE based calculations and here, the confinement relies on a 
A^.Aj ansatz, which does not respect local colour gauge invariance. Despite this defect, the 
first estimates have been done under this assumption, for simplicity. For baryons, the ansatz 
A^.Aj is a good approximation and it remains to be seen how good it is for pentaquarks. 

But, at any level of accuracy, it is expected that the upper bound obtained for pen- 
taquarks of negative parity, as those studied in Ref. [[UJ, would remain higher than the 
upper bound obtained for positive parity pentaquarks as introduced here. In an OGE based 
model, the situation will be the other way round. It remains to future experiments to 
disentangle between GBE and OGE models, through the search of pentaquarks. 
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TABLES 



TABLE I. Masses of hadrons required to calculate the threshold energy Et = Eb aryon +E meson . 
The experimental mass for mesons represents the average M = (M + 3M*) /4, when both the 
pseudoscalar M and the vector M* mesons masses are available |^0) . 



Hadron 


variational 


Mass (GeV) 


experiment 


N 


0.970 




0.939 


S 


1.235 




1.192 


D 


2.001 




1.973 


D s 


2.087 




2.076 


B 


5.302 




5.313 


B s 


5.379 




5.375 


Is 


175.711 
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TABLE II. Lowest positive parity pentaquarks of total angular momentum J= 1/2,3/2 and 
symmetry structure given by Eq. (2). Column 1 gives the flavour content, column 2 the isospin 
I, columns 3 and 4 the optimal variational parameters associated to the wave functions (15)-(18), 
column 5 gives AE = E — Et where E is the upper bound for the expectation value of the 
Hamiltonian (6)-(10) and Et = Et, aryon + E meson (see Table 1) for the threshold given in the last 
column. 

Pentaquark I Variational parameters (GeV 2 ) AE Threshold 

a b (MeV) 

uuddc 0.110 0.040 -75.6 N + D 

uuddb 0.110 0.053 -95.6 N + B 

uuddt 0.110 0.061 -102.8 N + ts 

uudsc 1/2 0.101 0.041 104.7 N + D S 



s 



uudsb 1/2 0.102 0.054 93.4 N + B 

uudst 1/2 0.102 0.063 86.9 N + ts 

udssc 1 0.092 0.041 81.0 E + D S 

udssb 1 0.092 0.055 69.3 S + S s 

udssi 1 0.093 0.064 62.6 S + is 



13 



